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ABSTRACT 

High-z Type la supernovae are expected to be gravitationally lensed by the foreground distribution of large- 
scale structure. The resulting magnification of supernovae is statistically measurable, and the angular correla- 
tion of the magnification pattern directly probes the integrated mass density along the line of sight. Measure- 
ments of cosmic magnification of supernovae therefore complements galaxy shear measurements in providing 
a direct measure of clustering of the dark matter. As the number of supernovae is typically much smaller than 
the number of sheared galaxies, the two-point correlation function of lensed Type la supernovae suffers from 
significantly increased shot noise. Neverthless, we find that the magnification map of a large sample of super- 
novae, such as that expected from next generation dedicated searches, will be easily measurable and provide 
an important cosmological tool. For example, a search over 20 sq. deg. over five years leading to a sample of 
^10,000 supernovae would measure the angular power spectrum of cosmic magnification with a cumulative 
signal-to-noise ratio of ^20. This detection can be further improved once the supernova distance measurements 
are cross-correlated with measurements of the foreground galaxy distribution. The magnification maps made 
using supernovae can be used for important cross-checks with traditional lensing shear statistics obtained in 
the same fields, as well as help to control systematics. We discuss two applications of supernova magnifica- 
tion maps: the breaking of the mass-sheet degeneracy when estimating masses of shear-detected clusters, and 
constraining the second-order corrections to weak lensing observables. 



Subject headings: cosmology: observations 
gravitational lensing 



cosmology: theory — galaxies: fundamental parameters 



1. INTRODUCTION 

Type la supernovae (SNe) are by now firmly established 
as powerful probes of the expansion history of the uni verse 
(Barri s et aljl2004t iKnop et alJl2003t iRiess et al] 120041) . In 
particular, the luminosity distance measurements from SNe 
provide a direct probe of dar k energy in the univ erse and its 
temporal behavior (e.g., Huterer & Coorav 2005 and refer- 
ences therein). Numerous current and future SNa la surveys 
are being planned or performed, and this community-wide ef- 
fort is expected to reach its apex with the NASA/DOE Joint 
Dark Energy Mission (JDEM). 

While SNe are very good (~ 10% errors in flux) standard 
candles, the inferred luminosity of a given supernova is af- 
fected by cosmic magnification due to gravitational lensing 
from the mass distribution of the large scale structure along 
the fine of sight between the supe rnova and the observer 
JFriemanll 19971 iHolz & Waldlll998l) . This is a fundamental 
limitation to the utility of standard candles, and SNe at high 
redshift (z > 1) are especially prone to fluctuations of their 
flux due to lensing. There was a recent claim of e vidence 
for weak lensing of SNe from the Riess et al. (2004) sample 
(Wang 2005); however, this claim remains unconfirmed as the 
correlation with foreground galaxie s that would be expecte d 
from lensing has not been observed (Menard & Dalai 2005). 

Weak lensing biases the luminosity measurement from each 
SN and thus introduces a systematic error in the extraction 
of cosmological parameters. With large number of SNe in 
each redshift interv al at high z this systematic can be essen - 
tially averaged out (Dal ai et al.ll2003t IHolz & Linderl l2005). 
though the full lensing covariance must be taken into account 
for accurate cosmological parameter estimates (Cooray et al. 
2005). While lensing has mostly been considered as a nui- 



sance, planned large area SN surveys provide an opportunity 
to treat lensing magnification on SNe as a signal, and extract 
information about the underlying dark matter distribution. In 
practice, by comparing the SN Hubble diagram averaged over 
all directions with individual SN luminosity distance mea- 
surements, one can map out anisotropy in the SN Hubble di- 
agram. This anisotropy will trace cosmic magnification, or 
in the weak gravitational lensing limit, it will be linearly pro- 
portional to the convergence and hence to fluctuations in the 
distribution of the foreground large scale structure. 

Previous studies have considered potential applications of 
cosmic converg ence as a p robe of the large-scale dark mat- 
ter distribution (Jain 2002). Cosmic magnification has al- 
ready been detected via cross-correlation between fluctua- 
tions in background source counts, such as quasars or X- 
ray sou rces, and a sample of low -redshift foreground galax- 
ies (see Bartel mann & SchneiderJ ( 120011) for a review). Past 
detections have often been affected by systematic uncertain- 
ties, but the powerful Sloan Digital Sky Survey (SDSS) has 
rece ntly ma de the first reliable detection of cosmic magnifica- 
tion (IScranton et alJl2005l) . Nevertheless, even this measure- 
ment using several hundred thousand quasars and upwards of 
ten million foreground galaxies was at a relatively modest 8cr 
level, illustrating the intrinsic difficulty in extracting the cos- 
mic magnification. Some proposals for future detections in- 
volve the use of 21 cm background a nisotropics of the hy dro- 
gen distribution prior to reionization ( Zhang & Pen 2003), but 
these studies are experimentally challenging and affected by 
large theoretical uncertainties in the amplit ude of the ex pected 
signal and its modification due to lensing dCooravl2004l) . 

In this Letter we propose mapping cosmic magnification 
with a sample of Type la SNe. We compute the predicted an- 
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FIG. 1. — Angular power spectrum of cosmic magnification (bottom 
long-dashed curve), cross-correlation between magnification and foreground 
galaxies (middle solid curve), and foreground galaxy clustering (top short- 
dashed curve). Fractional error in the luminosity of each SN has been as- 
sumed to be 0.1, and the error boxes account for both the sample (cosmic) 
variance due to the limited survey area and the presence of shot noise due to 
the finite number of SNe. We assume 10,000 SNe obtained over an area of 10 
sq. deg. with a uniform distribution in redshift between 0.1 and 1.7. The thin 
dashed line shows the shot noise in the magnification power spectrum. The 
cross-correlation makes use of the sample of SNe at z > 0.7, while the fore- 
ground galaxy sample is from Scranton et al. 1 2005). The thin dot-dashed line 
shows the second-order correction to magnification following Men ard et alj 
(2003) and using the halo model to describe the density field bispectrum. 



gular power spectrum of lensing magnification and estimate 
how accurately it can be measured, as well as how this mea- 
surement can be improved through cross-correlating super- 
nova distances with the foreground distribution of galaxies. 
We envision several important applications of this technique. 

The Letter is organized as follows: In §[2] we describe the 
weak lensing of SNe and the extent to which it can be mea- 
sured from future SN searches. In §[3] we discuss our results 
and comment on specific applications of SN magnification 
maps. We adopt the current concordance cosmological model 
with a Hubble constant of h = 0.7, matter density Q m = 0.3, 
cosmological constant A = 0.7, and normalization of the mat- 
ter power spectrum as = 0.85. 

2. WEAK LENSING OF SUPERNOVAE 

We begin by summarizing the effect of lensing magni- 
fication on SNe. While we concentrate on SNe, which 
have been firmly established as important and reliable cos- 
mological prob es, our study a pplies to any standard can- 
dle (e.g. lHolz & Hughes! J2005I) V Luminosity of a given su- 
pernova at a redshift z and located in the direction n, L(z,n), 
is affected by weak lensing magnification so that L(z,n) = 
n)L, where n) is the weak lensing induced magnifica- 
tion in the direction n and at redshift z, and L is the true lumi- 
nosity of the supernova. Note that /i can take values between 
the empty-beam value and infinity; the probability distribu- 
tion function of /i, P(/j,), has been extensively studied both 
analy tically and numerically (e.g., iHolzl 11998 : Wang et al. 
2002). Since (/i) = 1, one can av erage over large samples to 
determ ine the mean luminosity L ( Wang 20001 iHolz & Linderl 
120051) . One can now consider spatial fluctuations in the lumi- 
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which traces fluctuations in the cosmic magnification fi. In 
the weak lensing limit (/i, n <C 1) we have 

M = [(l- K ) 2 -|7| 2 r 1 «l+2« + 3« 2 +| 7 | 2 +..., (2) 



where k is the lensing convergence and \j\ = \Jjf + ^2 ls 
the total shear. To first order in the convergence, 6l(z,ti) 
traces spatial fluctuations of 2k, though high er order cor- 
rections may be important (Menar cfet 31.1120031) . Traditional 
weak lensing involves measurement of statistics of the shear, 
7,, as this leads to a distortion of background galaxy shapes 
(Bartelmann & Schneider 2001). Magnification, on the other 
hand, changes images sizes, but suffers from the problem that 
the true size of cosmological objects is highly uncertain. Fluc- 
tuations in the luminosity of standard candles provides a reli- 
able way to probe cosmic magnification. 

The angular power spectrum of magnification fluctuations 
is, assuming statistical isotropy 

(MfaMf'm') = CfSu'Smm' ■ (3) 

where fi£ m are the multipole moments of the magnification. 
Using the Limber approximation, the angular powe r spectrum 
can be written as (Kaiser 1998; Coorav et al. 2001) 
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l c 2 a(r) d A (r>) (4) 
where r is the comoving distance, d A is the angular diame- 
ter distance and n(r) is the radial distribution of SNe normal- 
ized so that J drn(r) = 1. Pd m (k,r) is the three-dimensional 
power spectrum of dark matter evaluated at the distance r; 
we calculate it using t he halo model of the large-scale struc- 
ture mass distribution (Coorav & Sheth 2002). The next or- 
der correction term, (K^ m K£> m >), is easily related to the con- 
vergen ce bispectrum and we calculate it using the halo model 
JCooravetalJ2001l) . 

In addition to a measurement of the projected angular 
power spectrum of cosmic magnification, one can also cross- 
correlate the magnification with the foreground galaxy distri- 
bution. The idea here is that the dark matter distribution that 
causes the magnification pattern <5/i(z,n) is traced by galax- 
ies and, therefore, 5^ and the normalized galaxy overdensity, 
5gai, are correlated. The projected cross-correlation between 
the two fields is described by the angular power spectrum 
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where n gd \(r) is the normalized radial distribution of fore- 
ground galaxies. Since 6i = 6^, the cross-correlation is 
independent of the power-law slope of the source number 
counts, unlike in the case of traditional galaxy-quasar cross- 
correlation measurements (Scranton et al. 2005). 

To estimate how well these angular power spectra can be 
measured with upcoming surveys, we compute the cumulative 
signal-to-noise ratio for detection 
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FIG. 2. — Left panel: Signal-to-noise ratio for the detection of the magnification power spectrum (dashed curves) and the magnification-galaxy cross power 
spectrum (solid curves) as a function of the smallest scale (maximum multipole) probed by the survey. We show cases of 3,000 SNe (cyan/light curves) and 
10,000 SNe (black/dark curves) collected over 20 sq. deg. (/^y m 0.0005). Note that the cross power spectrum can be detected with better signal-to-noise than 
the magnification auto-correlation because the shot noise in the foreground galaxy population is much smaller than that in the SNe number density. Right panel: 
Signal-to-noise ratio of the magnification power spectrum (black/dark curves) and the galaxy-magnification cross-power (orange/light curves) as a function of 
the fraction of the sky covered, /sky, and assuming 10,000 SNe (solid) and 100,000 SNe (dashed). Vertical lines show the minimal /sj™ in for a given number of 
SNe, which is given by the SN rate over the survey area and assuming an observing time of five years. 



where the index i references either the magnification power 
spectrum or the magnification-galaxy cross power spectrum. 
The error in the magnification power spectrum is given by 
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where A^sn is the surface density of SNe (number per stera- 
dian), is the uncertainty in the 5„ measurement from each 
supernova, / s k y is the fraction of sky covered by the survey, 
and A£ is the binning width in multipole space. For the SN 
luminosity-galaxy count cross correlation, the error is 
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where Cf ~ ga is the angular clustering power spectrum of 
foreground galaxies, and N ga \ is their surface density. 

For definitiveness we assume a magnification measurement 
error for each SN of = 0. 1 1 ; while smaller than errors in 
current SN surveys, this is expected to be achievable in the 
near future. For simplicity we consider A^sn SNe uniformly 
distributed in 0.1 < z < 1.7, which roughly approximates the 
distribution expected from SNAP ( Aldering et al. 2004). 

3. RESULTS AND DISCUSSION 

Figure^shows the angular power spectra of cosmic magni- 
fication of SNe, the cross power spectrum between the mag- 
nification anisotropies and the foreground galaxy distribution, 
and the galaxy angular power spectrum. Here we are partic- 
ularly interested in the and /i-gal spectra, for which we 

1 The magnification error is equal to the relative error in measuring lumi- 
nosity, and roughly equal to twice the relative luminosity distance error. 



also show error bars. We have assumed th e same foreground 
galaxy sample as in Scranton et al. (2005), with the redshift 
distribution of the form dn/dz ~ z L3 exp[-(z/0.26) 2 ' 17 ] and a 
number density of usable galaxies of 3 arcmin -2 . To avoid 
overlap of SNe with the galaxy sample, in the case of cross- 
correlation we only consider a SN subsample with redshifts 
greater than 0.7. In Figure^we also show the second order 
correction to magnification correspond ing to the k 2 ter m in 
Eq. (0; we calculate this following Menard et al. (20021), but 
using the halo model description for the density field bispec- 
trum rather than a numerical fit to simulations. 

In Figure |2 we show the total signal-to-noise ratio for the 
detection of the angular power spectrum of SN magnifica- 
tions. The left panel shows the S/N as a function of the small- 
est scale (maximum multipole) probed by the survey. We 
show cases of 3,000 and 10,000 SNe observed over 20 sq. deg. 
(i.e. / s ky ~ 0.0005). Note that, while the signal-to-noise ratio 
for the magnification power spectrum is usually around 20 or 
below, the cross power spectrum can be detected with consid- 
erably better significance due to the much smaller shot noise 
in the foreground galaxy population. Compare this to the cur- 
rent state-of-the-art in mapping the cosmic magnification: the 
SDSS catalog of quasars and galaxies has been used to detect 
the cosmic magnification at the 8<r level, and this detection 
is not expected to improve significantly given the already im- 
pressive statistics. In the future, SNe will provide the best 
opportunity to extend and improve the mapping of the cosmic 
magnification. The right panel of Figure [2] shows the signal- 
to-noise ratio, except now as a function of the sky coverage of 
the survey, f s ^, as we hold the total number of observed SNe 
fixed. This allows us to optimize the magnification measure- 
ment for a given amount of telescope time. Very small fay 
leads to large cosmic variance, while large / s k y decreases the 
surface density of SNe (since both their number and the ob- 
servation time are held fixed) and therefore increases the shot 
noise. In addition, the upper limit on the rate of SNe trans- 
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lates in to a minimum fi kv- Measurements of the actual rate 
of SNe flP ain et all2002) combined with theoretical estimates 
toda&Totani 2005) suggest that a year-long survey should 
find up to ~ 10 3 SNe per square degree, and this limit, as- 
suming a five-year survey, is shown as a lower limit on f^ y in 
the right panel of Figure|2] 

The surface density of galaxies, estimated to be around 
10 9 sr" 1 down to 27th magnitude (ISmail et alJll995l) . is far 
larger than the surface density of SNe (which is of order 10 6 
sr" 1 over a year-long integration). Nevertheless, cosmologi- 
cal methods that use galaxies to probe the formation of struc- 
ture in the universe, chiefly through weak gravitational lens- 
ing, are subject to systematic errors that range from theoreti- 
cal uncertainties to a variety of measurement systematics (see 
iHuterer et alJl20 05 and references therein). Supernova mea- 
surements of the magnification can be extremely valuable in 
helping control these systematics and break certain degenera- 
cies. 

For example, in order to establish the masses of shear- 
selected galaxy clusters, one can reconstruct the convergence 
k from the measured shear maps, 71 ?,(n) . While well-k nown 
techniques exist for this purpose (e.g. Kaiser & Sauires 1993), 
the reconstruction is insensitive to a constant mass sheet (or 
spatial ly uniform convergence). This ma ss-sheet degeneracy 
(e.g. lFalco et al.ll985llBradac et all2004l) can be broken with 
direct convergence measurements via magnification, and SNe 
are ideal candi dates for this purpose . Updating the calcula- 
tions in lKolatt & BartelmanrJ i ll 9981) . we find that up to 2% 
of the SNe are magnified by foreground clusters at a factor 
greater than 1.3 (a 3a or better detection). A survey covering 
~ 20 sq. deg. with ~ 10 4 SNe, combined with a shear analy- 
sis, can provide mass (enclosed out to the impact radius of the 
background supernova) of ~ 100 clusters to better than 10%. 

This approach can also test the consistency between shear 
measurements from galaxy shapes and convergence from SNe 
luminosity anisotropics. One can construct E- and B-modes 
of shear and, in the weak lensing limit, Cf = Cf and = 0. 
Departures from these relations are expected from both phys- 
ical and theoretical systematic uncertainties. For example, in- 
trinsic correlations between galaxies may pro duce add itional 
but u nequal contribution to E and B-modes l lHeavens et al.l 
2000). Moreover, there will exist contributions from higher- 
order effects due to slight departures from the weak lensing 
limit (see Eq. 0; Mena rdet alJl2003l) or higher order cor- 



rections to lensing that induce a rotational compon ent vi a 
coupling between two or more lenses dCoorav & Hull2002T) . 
or due to the presenc e of a gravitational wave background 
(Dodelson et al. 2003). The power of this consistency test is 
limited by the size of the higher order corrections. 

To quantify the detectability of the difference between the 
shear and convergence power spectra, we assume for a mo- 
ment that this difference is given by the second-order term 
Cg ~~ K plotted as the dot-dashed line in Fig. ^ We calculate 
the signal-to-noise in measuring the quantity Cf_= \C£ — CJ\ 
following the procedure similar to that in Eq. (|6j. Since fu- 
ture weak lensing shear surveys will have much smaller shot 
noise (7 r 2 ms /n ~ 10 -11 ) than the corresponding magnification 
power measurements (~ 10~ 9 ), the former source of noise can 
be ignored in the calculation. Using the SNe surface density 
of 10 3 deg" 2 yr" 1 , we find that this difference between the 
power spectra can be detected with a signal-to-noise ratio of 
~ 10(20 deg 2 /A)" 1,/2 , where A is the total survey area. Al- 
ternatively, if we assume that the fiducial difference between 
the power spectra has a shot-noise power spectrum (i.e. flat 
in £) with = A, the minimum detectable amplitude (with a 
signal-to-noise ratio of unity) is A w 5 x 10~ 7 (20 deg 2 /A) 1 / 2 . 
Consequently, corrections to the shear signal that are due to 
intrinsic correlations may be detectable (Jing 2002|). 

Magnification statistics from SNe also provide informa- 
tion on the cosmological parameters. This can be estimated 
in similar f ashion to the case of conventional w eak lensing 
of galaxies (Hu & Tegmark 1999; Huterer 2002), since both 
techniques probe the matter power spectrum and geometrical 
distance factors (see Eq. 0}. With the magnification power 
spectrum detected at a signal-to-noise ratio of 10 (100), one 
linear combination of parameters, typically with large weights 
in the fl m and directions, can be constrained to 10% (1%). 
While conventional weak lensing of galaxies can provide a 
more accurate overall determination of cosmological param- 
eters, the strength of the proposed method is that it combines 
lensing shear and magnification information in the same field, 
thereby providing a number of cross checks on systematics. 
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